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Abstract. The nearby galaxy M33 was observed by the 
imaging X-ray instruments on-board BeppoSAX. Two ob- 
servations at different phases of the 105.9 day intensity 
cycle of the luminous central source X-8 failed to reveal 
the expected modulation, suggesting that it is probably 
transitory. Similar behavior has been observed from sev- 
eral X-ray binary sources. This strengthens somewhat the 
idea that M33 X-8 is a black hole accreting from a binary 
companion. The 0.2-10 keV spectrum of M33 X-8 can best 
be modeled by an absorbed power-law with a photon in- 
dex, a, of 1.89 ± Q;yg and a disk-blackbody with a temper- 
ature, kT, of 1.10 ± 0.05 keV and a projected inner-disk 
radius of 55.4 ± f y km. This spectral shape is in good 
agreement with earlier ASCA results. The 2-10 keV spec- 
tra of M33 X-4, X-5, X-7, X-9 and X-10 are aU consistent 
with power-law or bremsstrahlung models, whilst that of 
X-6 appears to be significantly more complex and may 
be reasonably well modeled with a disk-blackbody with 
kT = 1.7 ± 0.2 keV and a projected inner-disk radius of 
7 ± 2 km. The spectrum of M33 X-9 is rather hard with 
a — 1.3. Compared to earlier Einstein and ROSAT ob- 
servations, M33 X-7 and X-9 varied in intensity, M33 X-4, 
X-6, and X-10 may have varied and M33 X-5 remained 
constant. 



Key words: Galaxies: individual (M33) 
- X-ray: galaxies 



Galaxies: nuclei 



1. Introduction 

M33 (NGC 598) is one of the nearest galaxies (795 kpc, 
van den Bergh 1991) and is classified as a late-type Sc spi- 
ral. As well as containing several point X-ray sources, M33 
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is unique amongst the Local Group in having a very bright 
(~10^^ erg s~^) X-ray source (X-8) located at a position 
consistent (to within the positional error of the Einstein 
High Resolution Imager (HRI) of 3") wi th the optical nu- 
cleus of the galaxy (see Long et al. 1981; Markert & Rallis 



1983; Gottwald et al. 1987; Trinchieri et al. 198?; Dubus et 



al. 1997). This source dominates the X-ray fiux from the 
galaxy, being responsible for ~70% of the 0.15-4.5 keV 
emission. Another X-ray source in M33 (X-7) is an eclips- 
ing X-ray binary pulsar (Peres et al. 1989; Schulman et al. 
1993|, [1994| Larson & Schulman |1997|) with a 3.45 day or- 



bital period and a 0.31 s pulse period (Dubus et al. 1999). 
The X-ray sources in M33 have been studied by Einstein, 
EXOSAT, ROSAT, and ASCA, as well as by earher non- 
imaging instruments. With an exposure of 50 ks ROSAT 
Position Sensitive Proportional Counter (PSPC) observa- 
tions of M33 have detected 35 sources in the central 30' 
which are brighter than 10^^ erg s~^ (Long et al. 1996; 
see also Haberl & Pietsch 2000). From an examination of 



archival ROSAT data spanning over 6 years (a total of 
387 ks of HRI time and 62 ks of PSPC time), Dubus et al. 
( 1997 ) find that the X-ray intensity of M33 X-8 is almost 
certainly modulated with a period of 105.9 ±0.1 days. 

The nature of M33 X-8 is unclear. X-ray transients 
such as Aql X-1 and 4U-1630-47 show outbursts on 
timescales of hundreds of days, but are not periodic (e.g., 



Kuulkers et al. 1997) while a number of persistent X-ray 



sources exhibit long term periodicities in their X-ray in- 
tensities. Extended observations using the All-Sky Moni- 
tor (ASM) on the Rossi X-ray Timing Explorer (R-XTE) 
have revealed changes in some long term periodicities (e.g., 
Paul et al. |2000[ ); Wojdowski et al. |1998| Kong et al. 
1998). The detection of a 106 day periodicity over ~20 



cycles implies a regularity and duty cycle incompatible 
with typical soft X-ray transients. M33 X-8 may instead 
be a moderate mass black hole. Until recently, the M33 nu- 
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Table 1. BeppoSAX M33 observing log. Phase, $, is determined using the 105.9 day ephemeris of Dubus et al. (1997) 
where the M33 X-8 maximum is expected at phase 0.0 and minimum at ~0.7 



Num. 


Observation 


Exposure 


M33 X-8 count rate 






Start 


End 


LEGS 


MECS 


LEGS 


MEGS 






(year mn dy hr:mn) 


(mn dy hr:mn) 


(ks) 


(ks) 


(0.1-2 keV; s"^) 


(1.8-10 keV; s"') 




1 


1998 Aug 24 13:05 


Aug 25 22:26 


32.3 


59.9 


0.068 ± 0.002 


0.157 ±0.002 


0.72 


2 


1999 Jan 06 20:16 


Jan 07 18:38 


0.9 


1.1 


0.045 ± 0.008 


0.131 ± 0.011 


0.98 


3 


2000 Jul 09 17:56 


Jul 13 04:33 


7.3 


37.4 


0.065 ± 0.004 


0.141 ±0.002 


0.20 



clear source was frequently considered (as are many of the 
bright sources in normal galaxies that exceed the Edding- 
ton limit of ~10^^ erg s~^ for a IMq compact object) to 
be a "mini-" (or perhaps "micro-" ) active galactic nucleus 
(AGN). However, the X-ray spectrum of X-8 appears to 
be much steeper than a classical AGN power-law and can 
be well fit by a multi-color dis k blac kbody model plus a 
harder power-law (Takano et al. 1994). This spectral shape 
is characteristic of the galactic black hole candidates which 
is often ultra-soft at low-energies and ultra-hard at high- 
energies. The ultra-soft component may originate from an 
optically thick accretion disk, while the ultra-hard com- 
ponent may result from comptonization of soft photons 



by hot electrons (Sunyaev & Triimper 1979). The accre- 
tion disk model is parameterized by only two quantities, 
the projected inner disk radius, rin(cosi)° '^, and the tem- 
perature at this inner radius, kTjn (Mitsuda et al. |1984| ; 
Makishima et al. |1986| ). Using ASCA, Takano et al. (^99§ 
obtain a value for rin(cosi)'''^ of ^^50 km, typical of black 
hole candidates in their high state. At the distance of M33, 
this corresponds to a black hole mass of '-^10 M©. 

In this paper we report on BeppoSAX observations 
of M33 primarily designed to investigate any spectral 
changes associated with the 105.9 day cycle of M33 X- 
8 (Sect. H) as well as the spectral and temporal properties 
of the luminous X-ray sources in M33 (Sect. 0). 



2. Observations 

Data from the coaligned Low-Energy Concentrator Spec- 
(LEGS; 



trometer (LEGS; 0.1-10 keV; Parmar et al. |1997D , 
Medium-Energy Goncentrator Spectrometer (MEGS; 1.8- 
10 keV; Boella et al. 1997), and the Phoswich Detection 



System (PDS; 15-300 keV; Frontera et al. |1997D on-board 
BeppoSAX are presented. The MEGS consists of two iden- 
tical grazing incidence telescopes with imaging gas scin- 
tillation proportional counters in their focal planes. The 
LEGS uses an identical concentrator system as the MEGS, 
but utilizes an ultra-thin entrance window and a drift- 
less configuration to extend the low-energy response to 
0.1 keV. The fields of view (FOV) of the LEGS and MEGS 
are circular with diameters of 37' and 56', respectively. 
The MEGS has a 2' thick circular window support struc- 
ture which is centered 10' from the FOV center. In addi- 
tion, four 2' thick radial spokes extend outwards from the 



circular support. The non-imaging PDS consists of four 
independent units arranged in pairs each having a sepa- 
rate collimator. Each collimator was alternatively rocked 
on- and 210' off-source every 96 s during the observations. 

M33 was observed 3 times by BeppoSAX (see Table p. 
The observing times were chosen in order to study M33 X- 
8 near the maximum and minimum of the 105.9 day X-ray 
cycle. The maxima occur at MJD = 48750.3 + n 105.9, 
and the minima approximately 0.7 in phase later (Dubus 
et al. 1997). The 1998 August observation was close to the 
105.9 day minimum ($ = 0.72), and the 1999 January ob- 
servation close to the maximum ($ = 0.98). However, this 
observation was curtailed due to technical problems and 
repeated in 2000 July at $ = 0.20, still close to the ex- 
pected maximum of the 105.9 day cycle. Good data were 
selected from intervals when the elevation angle above the 
Earth's limb was >5° and when the instrument configura- 
tions were nominal, using the SAXDAS 2.0.0 data analysis 
package. 

Fig. |l| shows the MEGS 1.8-10 keV image obtained by 
summing data from all 3 observations. The total exposure 
is 98.3 ks. For clarity, the data have been rebinned by a 
factor 2 into 16" pixels. No corrections for vignetting or 
obscuration by the MEGS window support structure (see 
above) have been applied. As well as the central M33 X-8 
source which dominates the overall X-ray flux, M33 X-4, 
X-5, X-6, X-7, X-9 and X-10 are clearly detected. The lu- 
minous source M33 X-3 is strongly obscured by the MEGS 
strongback and was not detected. In addition, the source 
identified here as M33 X-9 is resolved into 3 sources in 



ROSAT and iJmsiem observations (Long et al. 1996; Peres 
et al. |1989| ). 



3. M33 X-8 analysis 

LEGS and MEGS data were extracted centered on the 
position of M33 X-8 using radii of 6' and 4', respectively. 
This LEGS extraction radius is smaller than normally used 
in order to minimize the contribution from nearby sources. 
The cosmic and internal instrument backgrounds are ~5% 
of the source flux, and so background subtraction is not 
critical and was performed using source free regions of 
sky. Table |l| gives the 0.1-2.0 keV (LEGS) and 2.0-10 keV 
(MEGS) count rates for each observation. 
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Fig. 1. MECS 1.8-10.0 keV image of M33 obtained by summing data from all 3 observations. M33 X-8 is the central 
bright source. The data have been rebinned by a factor 2 into 16" pixels. The location of the sources discussed in Sect. ^ 
are indicated. The "cut-outs" at the top and bottom of the image are due to the removal of instrument calibration 
source events 



Table 2. MECS positions (J2000) of t he fai nt point sources (~1' uncertainty radius), compared with the ROSAT 
PSPC positions reported in Long et al. ( 1996 ). Source X-9 is unreso lved i n the MECS observation s, but resolved into 
three source (a,b,c) with ROSAT and Einstei n. For X-9 (Long et al. 1981 ) and X-9(c) (Peres et al. 1989) the Einstein 
co ordina tes are given. B74=Boulesteix e t al. 
al. 



1996; CS82a=Christian & Schommer 1982a; CS82b=Christian & Schommer 



1974| catalog of M33 H ii regions; L81— L ong et al. 

p82b| 



1981; L96=Long et 



Source 


R.A. (MECS) 


Dec. (MECS) 


R.A. (PSPC) 


Dec. (PSPC) 


Notes 




(h m s) 


(° ' ") 


(h m s) 


(° ' ") 




X-4 


01 33 17.5 


+30 54 02.0 


01 33 14.8 


+30 53 28.5 


SNR 136+396 (L96); H ii IC133 (L81) 


X-5 


01 33 25.9 


+30 44 33.8 


01 33 24.3 


+30 44 04.7 




X-6 


01 33 27.7 


+30 28 12.5 


01 33 27.5 


+30 27 22.7 


Globular cluster? (L81; CS82a) 


X-7 


01 33 31.8 


+30 33 36.4 


01 33 33.8 


+30 32 12.8 


X-ray eclipsing binary; H ll B208 (B74; L81) 


X-9 


01 34 26.0 


-1-30 55 59.7 


01 34 36.3 


-^30 55 00.4 


Einstein position (L81) 


X-9(a) 






01 34 38.7 


+30 55 10.3 




X-9(b) 






01 34 25.5 


+30 55 24.5 


Globular cluster? (n.3/4 in CS82b) 


X-9(c) 






01 34 31.5 


+30 56 20.5 


Elliptical galaxy at z=0.03? (CS82a) 


X-10 


01 34 50.8 


+30 29 17.3 


01 34 51.7 


+30 29 05.9 
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Table 3. Point sources detected in the MECS. C is the 
average MECS 1.8-10 keV background subtracted count 
rate. R is the extraction radius used for the MECS spec- 
tral analysis and 9 the source off-axis angle. Note that 
the MECS strongback obscures sources located between 
9 — 9-11'. V is the approximate vignetting correction to 
correct C to the on-axis case 



Source 


C 


R 


9 


V 




(10-3 s-^) 


(') 


(') 




X-4 


2.61 ±0.26 


4 


16 


2.4 


X-5 


3.40 ± 0.26 


2 


6 


1.3 


X-6 


6.15 ±0.29 


2 


12 


1.9 


X-7 


2.07 ±0.19 


2 


7 


1.3 


X-8 


151.8 ± 1.3 


4 





1.0 


X-9 


2.85 ±0.27 


4 


18 


2.7 


X-10 


8.95 ± 0.40 


4 


17 


2.5 



The M33 X-8 spectra were investigated by simultane- 
ously fitting data from the LECS and MECS instruments. 
All spectra were rebinned to oversample the full width 
half maximum of the energy resolution by a factor 3 and 
to have additionally a minimum of 20 counts per bin to 
allow use of the statistic. Data were selected in the en- 
ergy ranges 0.2-8.0 keV (LECS) and 1.8-10 keV (MECS) 
where the instrument responses are well determined and 
sufficient counts obtained. The photoelectric absorption 
cross sections of Morrison & McCammon (1983) are used 
throughout. Factors were included in the spectral fitting 
to allow for normalization uncertainties between the two 
instruments. 



3.1. 1998 August spectrum 

Initially, simple models were tried, including absorbed 
power-law, thermal bremsstrahlung and cutoff power- 
law (E"" exp — (Eco/kT)) models. The power-law and 
bremsstrahlung models give unacceptable fits with of 
380.8 and 149.8 for 101 degrees of freedom (dof), respec- 
tively. The cutoff power-law model fit is acceptable with a 
of 118.1 for 100 dof. Next, more complex models con- 
sisting of a power-law and a blackbody and a power-law 
and a disk-blackbody were fit. These both give acceptable 
fits with x^s of 106.0 and 108.9 for 99 dof, respectively. 
The best-fit inner disk radius is 55.2 ± 5 7 km for a dis- 
tance of 795 kpc and a disk inclination of 0°. The fit re- 
sults are summarized in Table ^. The absorption corrected 
1.0-10 keV source luminosity is 1.3 x lO'^^ erg s^^ for a 
distance of 795 kpc. 

3.2. 2000 July spectrum 

The spectra were created in the same way as for the 1998 
August observation and the same spectral models applied. 
The results are given in Table ^. The exposure times for 



this observation are shorter than for the 1988 August ob- 
servation, and the spectrum can be successfully modeled 
using either bremsstrahlung, cutoff power-law, power-law 
and blackbody or power-law and disk-blackbody models. 
The derived spectral parameters are consistent with those 
from the earlier observation with an inner disk radius 
of 50.0 ± 22 6 fo'^ tti^ same assumptions as above. 
The absorption corrected 1.0-10 keV source luminosity 
is 1.2 X 10'^^ erg s^^ for a distance of 795 kpc. 

3.3. Combined spectrum 

Since the spectral fit parameters derived from the fits to 
the 1998 August and 2000 July data are consistent, these 
data were combined to produce a "summed" spectrum. 
Data from the short (MECS exposure of 1.1 ks) observa- 
tion in 1999 January were also included. The same spectral 
models as above were fit and the results summarized in Ta- 
ble m As found by Takano et al. ( |1994[ ) , the best-fit model 
is the power-law and disk-blackbody combination with a 
kTin of 1.10±0.05 keV and an inner radius of 55.4 ± km 
for a distance of 795 kpc and a disk inclination angle of 0° . 
The mean, absorption corrected, 1.0-10 keV luminosity of 
M33 X-8 is 1.25 x 10^^ erg s^^ for a distance of 795 kpc. 
The 90% confidence upper limit to a narrow emission line 
at 6.5 keV is 100 eV. 

The galactic absorption in the direction of M33 X- 
8 is 5.6 x 10^° atom cm-^ (Dickey & Lockman |1990D. 



Using the power-law and disk-blackbody model, the best 
fit absorption of (1.2 ± Q;g) x 10^^ atom cm^^ suggests 
the presence of a small amount of additional absorption. 
The 3(7 upper-limit to any PDS counts in the 15-40 keV 
band is 0.06 s-^, consistent with the extrapolation of the 
best-fit power-law and disk-blackbody spectrum to this 
energy range which gives 0.02 count s"^ For a power-law 
spectrum with a = 2.1 the upper- limit corresponds to a 
flux of <0.4 mCrab. 



3.4. Time Variability 

Fig. I shows 2-10 keV MECS hghtcurves of M33 X-8 dur- 
ing 1998 August and 2000 July observations with a bin- 
ning time of 1024 s. The LECS and MECS data from these 
observations were searched for the presence of periodic 
signals. None were found. During the 1998 August obser- 
vation the 3cr limit to any 2-10 keV (MECS) root mean 
square variability is <7.6% when evaluated on a timescale 
of 3000 s. In the energy band 0.1-2.0 keV (LECS) the 
upper limit evaluated in the same way is <5.3%. For the 
2000 July observation the corresponding upper limits are 
<25% and <29%. 

The mean ROSAT HRI count rate is 0.20 s'^ and 
the modulation appears sinusoidal-like with an amplitude 
of ~20%. The folded 0.1-2.5 ke V ROS AT hghtcurve (see 
Fig. ^ and Fig. 3 of Dubus et al. 1997 ) has a mean count 
rate of 0.20 s^^ with a sinusoidal-like modulation super- 
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Table 4. M33 X-8 fit results. Nh is the absorption in units of 10^^ atom cm ^. rin(cosi)'''^ is in units of km for a 

distance of 795 kpc. a is the power-law photon index. 90% confidence limits are given when the per dof is <2 

Model Nh kT/kTin (keV) a Eco (keV) ^^(cosi)'''^ xVdof 
1998 August observation 

Power-law 0.58 ... 2.7 ... ... 380.8/101 

Bremsstrahlung 0.24 ± 0.03 2.93 ±0.10 ... ... ... 149.8/101 

Cutoff power-law 0.12 ±0.04 ... 0.63 ± 0.29 1.89 ±0.17 ... 118.1/100 

Power-law ± blackbody 0.23 ± 0.05 0.75 ± 0.03 2.44 ±0.15 ... ... 106.0/99 

Power-law ± disk-blackbody 0.15 ± 0.06 1.12 ±0.06 2.15 ± S;*!] ... 55.2 ± 108.9/99 

2000 July observation 

Power-law 0.55 ±0.10 ... 2.60 ± 0.07 ... ... 105.3/71 

Bremsstrahlung 0.21 ± 0.06 3.40 ± 0.20 ... ... ... 75.4/71 

Cutoff power-law 0.25 ±0.10 ... 1.49 ± 3.79 ± J;^^ ... 75.2/70 

Power-law ± blackbody 0.29 ±0.15 0.72 ± 2.33 ± 0.20 ... ... 76.9/69 

Power-law ± disk-blackbody 0.20 ± 0.20 1.05 ± o;?J 2.06 ± gi^g . . . 50.0 ± 74.0/69 

All observations 

Power-law 0.55 ... 2.7 ... ... 396.6/105 

Bremsstrahlung 0.22 ± 0.03 3.10 ±0.09 ... ... ... 141.9/105 

Cutoff power-law 0.135 ±0.05 ... 0.85 ±0.18 2.22 ± g:! ... 125.4/104 

Power-law ± blackbody 0.25 ± 0.05 0.75 ± 0.03 2.40 ±0.12 ... ... 122.2/103 

Power-law ± disk-blackbody 0.12 ± 1.10 ±0.05 1.89 ± arSJ ... 55.4 ± ?;? 115.7/103 



BeppoSAX LECS + MECS M33 X-8 BeppoSAX LECS + MECS M33 X-8 




0.2 0.5 12 5 10 0.2 0.5 12 5 10 

Energy (keV) Energy (keV) 

Fig. 2. The M33 X-8 LEGS and MEGS count spectrum and residuals when the power-law and disk-blackbody model 
is fit to the summed spectrum (left panel). The residuals are given in count s^^ keV^^. The right panel shows the 
photon spectrum with the contribution of the power-law indicated 



posed with an amplitude of ~20%. This implies that the 
count rates of the 1999 January and 2000 July observa- 
tions should be '^40% higher than that of 1998 August, 
should the 105.9 day modulation be present. However, the 
MEGS 2.0-10.0 keV count rate at $ = 0.72 (close to the 
expected minimum of the cycle) is higher than at the ex- 
pected maximum. In order to compare our BeppoSAX re- 



sult s for the two long observations with those of Dubus et 
al. ( 1997 ), we folded the best-fit BeppoSAX power- law and 
disk-blackbody models (see Table ^) through the HRI on- 
axis response matrix obtained from the HEASARG. This 
gives predicted HRI count rates of 0.221 ± 0.007 s~^ and 
0.195 ± 0.012 s^i for the 1998 August and 2000 July ob- 
servations, respectively. These points are shown together 
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Time (h) 




Fig. 3. MECS 2-10 keV lightcurves for M33 X-8 with a 
binning of 1024 s for the 1998 August (upper panel) and 
2000 July (lower panel) observations. No background has 
been subtracted 




105.9 day Phase 



Fig. 4 . The ROSAT HRI count rates of Dubus ct al. 
{ 09T\ ) folded over the 105.9 day period. The count rates 
corresponding to the two BeppoSAX measurements are 
indicated with filled circles 



with the HRI results of Dubus et al. ( |l997D on Fig. |. The 
derived count rates do not support the continued pres- 
ence of the 105.9 day modulation, although we cannot en- 
tirely exclude the possibility that the modulation is being 
masked by a higher than usual level of intrinsic variability. 

4. Point sources 

In order to investigate the properties of the other sources 
identified in Fig. MECS data from all 3 observations 
were used. Except for M33 X-6 LEGS data were not used 
since there were too few counts present. During the 1998 
August observation M33 X-4 was partially obs cured by a 
MECS strongback radial spoke (see Sect. L2) and these 



data were excluded from the analysis. For all the other 
sources no obscuration problems were present. Spectra 
were extracted centered on each source using the extrac- 
tion radii listed in Table ^ and 2' in the case of the LEGS 
M33 X-6 data. These radii were chosen to minimize con- 
tamination from nearby sources and take into account the 
position of the MECS window support strongback. The 
spectra were rebinned in the same way as in Sect. ^ and 
data selected in the energy ranges where sufficient counts 
were obtained. Since these sources are all faint by Bep- 
poSAX standards care was taken to reliably estimate the 
background. In order to properly take into account all 
the possible contaminating contributions, especially from 
the bright nuclear source M33 X-8, local backgrounds, 
each extracted from the same off-axis angle as the cor- 
responding source, but at the geometrically opposite posi- 
tion in the FOV, were used. This is also necessary because 
the MECS instrument background has a component due 
to Mn Ka calibration source X-rays which are detected 
within the FOV. This effect is most noticeable close to 
the calibration sources and in particular for M33 X-9. A 
comparison of these local backgrounds with the average 
properties of background taken from high latitude point- 
ings reveals that no strong evidence for diffuse emission 
is present in the MECS M33 data. A higher local back- 
ground (~20%) is present around M33 X-8, but this is 
consistent with being due to the X-ray mirror scattering 
wings. 

Instrument response files appropriate to the off-axis 
positions of each source were used, except for M33 X-5 
and X-7, which are close enough to the on-axis position 
to use the standard response files. 

4-1. Spectral fit results 

Absorbed power-law and bremsstrahlung models were first 
fit to each of the spectra and the results presented in Ta- 
ble ^. Fits were repeated with the absorption constrained 
to be greater than, and equal to, the galactic value for 
each source. With the exception of M33 X-6 (and pos- 
sibly M33 X-10) these models give reasonable fits to all 
the spectra and it is not possible to exclude either model. 
In the case of M33 X-6, the best-fit absorbed power-law 
and bremsstrahlung models give x^'s of 51.7 and 44.2 for 
28 dof, respectively. The best-fit to the M33 X-6 spec- 
trum is instead obtained with an absorbed disk-blackbody 
model with kTjn = 1.7 ± 0.2 keV and a normalization 
(rin(cosi)''-^) of 7±2 km. The 90% confidence upper limit 
to a narrow emission line at 6.5 keV from M33 X-6 is 
540 eV. In the case of M33 X-10, a broad feature appears 
to be present at keV. This can be modeled by the 
addition of a (cr = 200 eV) Gaussian line to the best-fit 
power-law to give a of 18.6 for 15 dof. The line energy 
is 4.2±g;2 keV and the equivalent width 440 eV. 

Table || includes estimates of the source luminosi- 
ties extrapolated (when necessary) to the 2-10 keV en- 
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Table 5. Spectral fit results to the M33 point sources. The absorption, Nh, is in units of 10^^ atom cm ^. Uncertainties 
arc quoted at 90% confidence. The 2-10 keV absorption corrected source luminosity, L, assumes a distance of 795 kpc 

Source Model Nh kT q L xVdof 

(keV) (erg s'^) 



X-4 


Power-law 


2.0 ± f:i 






3.4 ± 


2.8 
1.0 


f .5 X 10^* 


i3.7/f2 




Power- law 


0.056 fixed 






2.7 ± 


0.6 
0.4 


7.9 X fO^'' 


f5.0/i3 




Bremsstrahlung 


<4.9 


f.8 ± 


1.8 
0.9 






8.6 X iO^^ 


i3.2/f2 




Bremsstrahlung 


0.056 fixed 


2.3 ± 


1.4 
0.6 






7.0 X iO^'' 


f3.7/f3 



X-5 


Power-law 


<5.6 






3.2 ± 


1.2 
0.9 


3.4 X fO^'^ 


f0.7/i4 




Power- law 


0.056 fixed 






2.6 ± 


0.4 
0.4 


2.7 X iO^^ 


ff.6/f5 




Bremsstrahlung 


<2.9 


2.9 ± 


1.5 
1.2 






2.6 X iO^'^ 


9.7/f4 




Bremsstrahlung 


0.056 fixed 


2.9 ± 


1.7 
0.7 






2.5 X iO^'^ 


9.7/f3 



X-6 Power-law 


0.45 ± III 




2.1 ± 


0.2 
0.3 


7.6 X 10^'' 


51.7/28 


Power-law 


0.056 fixed 




1.8 ± 


0.1 
0.1 


7.5 X 10^'' 


56.9/29 


Bremsstrahlung 


<0.69 


5.6 ± li 






7.1 X 10^'^ 


44.2/28 


Bremsstrahlung 


0.056 fixed 


6.4 ± li 






7.1 X 10^^ 


45.0/29 


Blackbody 


<0.19 


f.OiO.lO 






7.1 X 10^^ 


41.4/28 


Disk-blackbody 


<0.29 


f.7±0.2 






6.8 X 10^^ 


37.2/28 



X-7 


Power-law 


<9.8 




2.9 ± 


1.7 
1.3 


2.8 X 10^'' 


6.2/9 




Power-law 


0.056 fixed 




1.7 ± 


0.6 
0.6 


1.6 X 10^'' 


7.3/10 




Bremsstrahlung 


<11 


3.7 ± 






1.8 X 10^'^ 


6.1/9 




Bremsstrahlung 


0.056 fixed 


>3.7 






1.4 X 10^'^ 


6.9/10 



X-9 


Power-law 


<3.9 




1.3 ± 


2.3 
0.8 


6.1 X 10^'' 


7.2/8 




Power-law 


0.056 fixed 




1.1 ± 


1.0 
0.6 


6.8 X 10^'' 


7.5/9 




Bremsstrahlung 


<4.2 


>2.6 






6.1 X 10^'' 


7.2/8 




Bremsstrahlung 


0.056 fixed 


>5.3 






5.7 X 10^'' 


7.2/9 



X-10 Power-law 


<3.3 






1.9 ± 


0.8 
0.4 


1.4 X 10^** 


31.2/18 


Power- law 


0.056 fixed 






1.7 ± 


0.2 
0.3 


1.4 X 10^** 


31.9/19 


Bremsstrahlung 


<2.0 


8.7 ± 


14 
4.6 






1.4 X 10^* 


31.1/18 


Bremsstrahlung 


0.056 fixed 


11 ± 


13 
4.6 






1.3 X 10^* 


31.2/19 



ergy range. Photon spectra, deconvolved using the best-fit 
power-law models are shown in Fig. I, except for M33 X-6 
where the best-fit disk-blackbody model is used. With the 
exception of M33 X-9 and X-10, all the sources appear to 
be relatively soft in the 2-10 keV energy range. 

4-. 2. Time variability 

In order to investigate if any of the point sources exhibit 
excess variability, a Kolmogorov-Smirnov test as imple- 
mented in the FTOOLs package lcstat was performed 
on each of the lightcurves presented in Figs. ^ and ^ The 
sources M33 X-4, X-5, and X-10 show some evidence for 
excess variability on time scales >1024 s with probabilities 
for being constant of 4.8%, 4.4%, and 5.4% for M33 X- 
4, (obs. 1), X-5 (obs. 3), and X-10 (obs. 1), respectively. 
The other 3 sources do not show evidence for variability 
during either long observation with probabilities for being 
constant of 32% (X-6, obs. 1), 61% (X-7, obs. 1, >2048 s), 
and 19% (X-9, obs. 3). For comparison, data from source 
free regions of the image were also extracted and processed 



in the same way. No evidence for any excess variability was 
found. 

A comparison between the 1998 August (Fig. ^), and 
2000 July (Fig. ^ observations, does not show evidence 
for variability, except for M33 X-4, which showed a fac- 
tor of '^2 difference in count rate between the two ob- 
servations. However, this change is almost certainly due 
to obscuration by one of the radial spokes of the MECS 
strongback during the 1998 August observation. The 2000 
July observation was made at a different spacecraft roll 
angle, resulting in the source being well clear of the strong- 
back. The BeppoSAX data are therefore consistent with 
the association of M33 X-4 with the s upern ova remnant 
SNR 136-^396 discussed in Long et al. (|l996|) . 

Table H gives a comparison of the BeppoSAX source in- 
tensities with those derived from previous ROSAT PSPC 



(Long et al. 1996) and Einstein IPC observations (Peres 



et al. 1989, For M33 X-4 only the data from the 2000 July 



observation is used. Unfortunately, due to the lack of good 
estimates of the absorbing column to most of the sources 
(see Table ||), it is difficult to reliably extrapolate the 
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Fig. 5. MECS photon spectra of the pomt sources detected in the M33 galaxy deconvolved using an absorbed power- 
law model except M33 X-6 where the LEGS spectrum is included and an absorbed disk-blackbody model is used. Data 
from all 3 observations are included, except for M33 X-4 where the 1998 August observation is excluded. The best-fit 
parameters are given in Table |^. The plots all have the same extrema to aid comparison 



BeppoSAX spectra to the PSPC (0.1-2.4 keV) and IPC 
(0.2-4.5 keV) energy ranges. Instead, we have adopted the 
spectral model assumed by Long et al. (1996) to convert 
ROSAT PSPC counts to fluxes. This consists of a power- 
law with a = 2 and absorption of 10^^ atom cm~^. We 
note that none of the spectral fits reported in Table |^ are 
in strong disagreement with such a model. We estimate 
that the uncertainty in the BeppoSAX fluxes is ~40%, 



primarily due to uncertainties in instrumental calibration 
at large off-axis angles. 

Examination of Table || shows that M33 X-4, X- 
6 {Einstein only), X-7, X-9 and X-10 (ROSAT only) 
show >40% variability with respect to both the Einstein 
and ROSAT measurements. The source M33 X-7 shows 
the strongest evidence for variability being a factor ~3 
fainter during the BeppoSAX observations than previ- 
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Fig. 6. Lightcurves for the 6 faint point sources during the 1998 August observation. No background has been sub- 
tracted 



ously observed. Since this source exhibits 0.33 day dura- 
tion echpses, we verified that the observations were not at 
the expected times of echpse using the ephemeris of Larson 
& Schulman ( |l997|) . Source X-9, not resolved in the MECS 
image, appears to be variable being a factor '^2.5 brighter 
during the BeppoSAX observation than previously, if the 
ROSAT fluxes of the three separate sources are summed 
toget her (the X-9(c) flux is not reported in Long et al. 



1996, although it is visible in their figure and appears to 
be as bright as the other two sources). Since the level of 
variability of M33 X-4, X-6, and X-10 is comparable with 
the (hard to quantify) systematic uncertainties, we regard 
the detection of variability for these sources as tentative. 

Since M33 X-4 may be associated with a SNR, 
it is interesting to compare BeppoSAX fluxes with 
those obtained using Einstein and ROSAT using 
a bremsstrahlung model. If Nh is again fixed at 
10^1 atom cm-^ we obtain ¥1=2 A keV for the MECS 
spectrum and 0.1-2.4 keV and 0.2-4.5 keV luminosities 
of 1.1x10^^ ergs s~^ and l.SxlO'^® ergs s~^, respectively. 
These values should be compared with 9.5x10"^'' ergs s^^ 
{Einstein IPC, 0.2-4.5 keV) and 8.3x103"^ ergs s~^ 
(ROSAT, 0.1-2.4 keV), derived using this same spectral 
shape. Thus, there is no strong evidence for flux variability 
from M33 X-4. 



Table 6. Comparison of the Einstein IPC and ROSAT 
PSPC, luminosities with those derived by extrapolating 
the BeppoSAX values to the Einstein (0.2-4.5 keV) and 
ROSAT (0.1-2.4 keV) energy ranges. See the text for de- 
tails of the assumed spectrum. The absorption corrected 
luminosity, L, assumes a distance of 795 kpc and is in units 
of 10^"^ erg s~^ 



Source 


L (0.1-2.4 keV) 


L (0.2-4.5 keV) 




BeppoSAX 


ROSAT 


BeppoSAX 


Einstein 


X-4 


18.9 


11.1 


18.6 


12.4 


X-5 


5.80 


5.72 


5.71 


6.84 


X-6 


14.0 


18.2 


8.57 


16.4 


X-7 


2.8 


12.3 


2.79 


9.00 


X-9 


8.6 




8.57 




X-9(a) 




1.7 




1.80 


X-9(b) 




1.9 




0.67 


X-9(c) 








0.66 


X-10 


29.0 


18.4 


23.6 


16.6 



5. Discussion 

The nearby galaxy M33 was observed by the imaging X- 
ray instruments on-board BeppoSAX. As well as the cen- 
tral M33 X-8 source which dominates the overall X-ray 
flux, M33 X-4, X-5, X-6, X-7, X-9 and X-10 were clearly 
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Fig. 7. Lightciirves for the 6 faint point sources during the 2000 July observation. No background has been subtracted 



detected. Observations close to the expected maximum 
and minimum of the 105.9 day intensity cycle of M33 X-8 
failed to reveal the expected modulation, suggesting that 
it is probably transitory. No obvious spectral differences 
between the observations are present. The 0.2-10 keV 
spectrum of M33 X-8 can best be modeled by an absorbed 
power-law with a — 1.89 ± g;|g and a disk-blackbody with 
kT = 1.10 ± 0.05 keV and a projected inner-disk radius 
of 55.4 ± 7 Y km. This spectral shape is in g ood a gree- 
ment with the ASCA results of Takano et al. (1994) and 
corresponds to a black hole mass of ~10 Mq. 

The nature of M33 X-8 is intriguing. X-ray tran- 
sients such as Aql X-1 and 4U 1630-47 show outbursts 
on timescales of hund reds of days, but are not periodic 
(e.g., Kuulkers et al. 



1997) 



A number of persistent X- 
ray sources exhibit long term periodicities in their X-ray 
intensities. Perhaps the best known of these are the pul- 
sars Her X-1 (35 days, Giacconi et al. |1973D, LMC X-4 



(30.5 days, L ang et al. |198lD and SMC X-1 (60 days, Wo- 
jdowski et al. 1998). In addition, the high-mass X-ray bi- 
nary SS 433 thought to contain a black hole exhibits a 
164 day photometric and spectroscopic modulation in the 
optical band (Margon et al. 1979). Two other black hole 



candidates, Cyg X-1 (300 days, Priedhorsky et al. 1983 
and LMC X-3, (198 or 99 days, Cowley et al. |1991|) also 



exhibit long term X-ray modulations. Thus, the presence 
of a 105.9 day X-ray intensity modulation from M33 X-8 is 



not particularly unusual if the source does indeed contain 
a '-lO M© black hole. 

Recently, extended observations using the ASM on R- 
XTE have revealed that in at least four sources (Cyg X-3, 

Paul ct al. ^000 ; 



LMC X-3, SMC X-1, and GX 354-0 

Wojdowski et al. 1998; Kong et al. 199^ ) the long term 



modulations have disappeared, or changed periods. Such 
unstable behavior is difficult to reconcile with the pre- 
ferred model for the modulations in the high-mass sources 
which is obscuration by a precessing accretion disk (see 
e.g., Priedhorsky & Holt 1987). Since the detection of 



changes in the properties of these long term variations 
requires very long observing baselines, it is possible that 
such unstable behavior is relatively common in accreting 
systems. Thus, the probable cessation of the 105.9 day 
periodicity reported here strengthens somewhat the idea 
that M33 X-8 consists of a black hole accreting from a 
binary companion. 

The 2-10 keV spectra of M33 X-4, X-5, X-7, X-9 and 
X-10 are all consistent with power-law or bremsstrahlung 
models, whilst that of X-6 appears to be significantly more 
complex, and may be reasonably well modeled with a disk- 
blackbody with kT = 1.7 ±0.2 keV and a projected inner- 
disk radius of 7 ± 2 km. This normalization is in good 
agreement with the average of a number of low-mass X-ray 
binary systems known to contain neutron stars (~7 km), 
but is smaller than the average (~25 km) for black hole 
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candidates (e.g., Tanaka & Lewin (1995). This suggests 
that M33 X-6 may be a neutron star accreting from a 
low-mass companion. The spectrum of M33 X-9 is quite 
hard with a = 1.3, similar to the hard/low state in galactic 
black hole candidates. A comparison with earlier Einstein 
and ROSAT results implies that M 33 X-7 and X-9 are 
variable on long timescales, M33 X-4, X-6 and X-10 may 
be variable, whilst X-5 does not show evidence for any such 
variability. We cannot exclude the association of M33 X-4 
with the su pernov a remnant SNR 136-1-396 discussed in 
Long et al. (|1996|) . The lack of variability, moderate lu- 
minosity (3 X 10'^^ erg s~^), and simple (consistent with 
a power-law or bremsstrahlung) spectrum may indicate 
that M33 X-5 is a young X-ray bright supernova rem- 
nant. M33 X-6 clearly has an unusual spectrum and is 
worthy of further study with more sensitive instruments. 
It is unsurprising that M33 X-7 shows strong long term 
variability since it is know n to b e an (eclipsing) X-ray bi- 
nary system (Peres et al. |989|). Whilst M33 X-9 was a 
factor ~2.5 brighter during the BeppoSAX observations 
than previously observed, it is unclear which (if any) of 
the 3 sources that apparently comprise X-9 had become 
significantly brighter. 
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